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ABSTRACT 

We study the effects of feedback from active galactic nuclei (AGN) on emission from molecular gas in 
galaxy mergers by combining hydrodynamic simulations which include black holes with a three-dimensional, 
non-local thermodynamic equilibrium (LTE) radiative transfer code. We find that molecular clouds entrained 
in AGN winds produce an extended CO morphology with significant off-nuclear emission, which may be 
detectable via contour mapping. Furthermore, kinematic signatures of these molecular outflows are visible in 
emission line profiles when the outflow has a large line of sight velocity. Our results can help interpret current 
and upcoming observations of luminous infrared galaxies, as well as provide a detailed test of subresolution 
prescriptions for supermassive black hole growth in galaxy-scale hydrodynamic simulations. 
Subject headings: cosmology: theory — galaxies: formation — active — interactions — ISM — line: formation 



1. INTRODUCTION 

The physical processes giving rise to the birth and sus- 
tained fueling of massive starbursts and active galactic nu- 
clei (AGN) have been of interest since their discovery. It is 
generally agreed that AGN are powered by accretion of gas 
onto supermassive black holes (BHs) in the centers of galax- 
ies (e.g. Lynden-Bell, 1969), but the fueling mechanism is 
less clear. Hydrodynamic simulations have shown that merg- 
ers can produce strong, galaxy-scale inflows owing to gravi- 
tational torques (Barnes & Hernquist 1991, 1996), triggering 
starbursts (Mihos & Hernquist, 1996). This suggests a cir- 
cumstantial link between starbursts and AGN activity. 

Observational evidence connecting the evolution of star- 
bursts into AGN through mergers is compelling. Some ultra- 
luminous infrared galaxies (ULIRGs, Lir >10 12 L©) exhibit 
spectral energy distributions (SEDs) characteristic of classi- 
cal starbursts, whereas others have SEDs more closely resem- 
bling optical quasars (e.g. Farrah et al. 2003). In seminal 
papers, Soifer et al. (1987), and Sanders et al. (1988a,b) used 
optical, near infrared and millimeter-wave observations to ad- 
vance a scenario in which starburst dominated ULIRGs served 
as antecedents of AGN. The specifics of a starburst- AGN con- 
nection, however, remain under heavy debate throughout the 
literature (e.g. Sanders & Mirabel, 1996). 

Recent numerical models by Di Matteo, Springel & Hern- 
quist (2005), Hopkins et al. (2005a-d), and Springel, Di Mat- 
teo & Hernquist (2005a) have provided a theoretical foun- 
dation for the link between starbursts and AGN. In particu- 
lar, by modeling the growth of (and feedback from) central 
black holes, they showed that gas-rich galaxy mergers are a 
viable candidate to serve as a precursor to the formation of 
quasars. Their simulations also show that feedback from ac- 
creting black holes in galaxies are relevant to a wide range of 
phenomena associated with the evolution of galaxies in merg- 
ers including characteristic X-ray emission patterns (Cox et 
al. 2006), observed quasar luminosity functions and lifetimes 
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(Hopkins et al. 2005c), the M B h - cr relation (Di Matteo et 
al 2005; Robertson et al. 2006a), the fundamental plane of 
ellipticals (Robertson et al 2006b), and the bimodal galaxy 
color distribution (Springel, Di Matteo & Hernquist, 2005b; 
Hopkins et al. 2005e). 

There has been a longstanding interest in better understand- 
ing the nature of the molecular interstellar medium (ISM) in 
mergers, as the molecular gas serves as fuel for the induced 
starburst activity and possibly for accreting BH(s). High reso- 
lution observations have identified massive concentrations of 
molecular gas in the nuclear regions of ULIRGs (Bryant & 
Scoville, 1999), as well as high excitation molecular gas in 
regions of massive starbursts (Iono et al. 2004; Wang et al. 
2004). 

Large surveys at submillimeter (sub-mm) and millime- 
ter wavelengths have shown that mergers in the local Uni- 
verse emit copious molecular line radiation both from diffuse 
molecular gas (e.g. Sanders, Scoville & Soifer, 1991, among 
others), and dense cloud cores (e.g. Gao & Solomon, 2004; 
Narayanan et al. 2005). Other work demonstrates that high- 
Z infrared luminous and sub-mm selected sources also con- 
tain significant amounts of molecular gas (Greve et al. 2005; 
Tacconi et al. 2006). Moreover, as evidenced by IR and X- 
ray studies, a large fraction of these galaxies at high-z con- 
tain AGN (e.g. Alexander et al. 2005, Polletta et a. 2006). 
However, despite the wealth of data on molecular emission in 
interacting galaxies, little is known about the impact of em- 
bedded AGN on this radiation. 

Here, we describe preliminary attempts to quantify the ob- 
servable effects of AGN feedback on molecular line emission 
from major galaxy mergers. We use hydrodynamic simula- 
tions of mergers with and without AGN, combined with a 
new 3D non-LTE radiative transfer code (Narayanan et al. 
2006a,b) to model CO emission. In this Letter, we describe 
results in which we find distinct signatures of AGN feedback 
on cold gas, and discuss some observational results that may 
be understood in this context. 

2. NUMERICAL SIMULATIONS 

Our hydrodynamic simulations were performed using 
the iV-body/smoothed particle hydrodynamics (SPH) code, 
GADGET-2 (Springel, 2005). This code uses a fully conser- 
vative formulation of SPH (Springel & Hernquist 2002), and 
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accounts for radiative cooling of the gas (Dave et al. 1999), 
a multi-phase description of the ISM which includes cold 
clouds in pressure equilibrium with hot, diffuse gas (e.g. Mc- 
Kee & Ostriker, 1977), and a prescription for star formation 
constrained by the Schmidt/Kennicutt laws (Kennicutt, 1998, 
Schmidt, 1959; see Springel & Hernquist, 2003). The black 
hole(s) in the simulation are realized through sink particles 
which accrete gas from the surrounding ISM such that 0.5% 
of the accreted mass energy onto the central black hole(s) is 
reinjected into the ISM as thermal energy (Di Matteo et al. 
2005, Springel et al. 2005a,b). 

The progenitor disk galaxies used in this work are simi- 
lar to the Milky Way, but with a higher gas fraction (50%). 
Thus, they are likely representative of high-z disk galaxies, 
and possible progenitors of present epoch ULIRGs. The 
methodology for constructing the model galaxies is given in 
Springel et al. (2005a). The progenitors utilized a softened 
equation of state (EOS) with softening parameter ^eos=0.25 
(Springel et al. 2005a) such that the mass-weighted ISM 
temperature is ^10 45 K. The galaxies had dark matter ha- 
los initialized to follow a Hernquist (1990) profile, and cir- 
cular velocity V2oo=160 km s" 1 . The virial properties of the 
halos are scaled to be appropriate for z = 2 (e.g. Robertson 
et al. 2006a), and follow the prescription given by Mo et al. 
(1998) for cosmological models of disk galaxies. The galax- 
ies were set on a parabolic orbit with the orientation of the 
spin axis of each disk specified by the standard spherical co- 
ordinates, 9 and </> (with 6>i=30°, </>i=60°, <9 2 =-30°, </> 2 =45°), 
and were initially separated by 140 kpc . We utilize 120,000 
dark matter particles, and 160,000 total disk particles, 50% 
of which represent gas, the rest serving as collisionless star 
particles. The gas, star, disk and dark matter particle masses 
were 3.9 x 10 5 , 1 .95 x 10 5 , 5.9 x 10 5 and 7.6 x 10 6 It'Mq each, 
respectively. The gravitational softening lengths were 100 
fr'pc for baryons, and 200 h _1 pc for dark matter particles. 
We have performed simulations with and without black holes. 
For the simulations with BHs, the initial mass of the BH par- 
ticle in each progenitor 10 5 h _1 M Q , and the peak accretion 
rate was ^0.5 M yr _1 . The final mass of the remnant's BH is 
-SxIO^-'Mq . 

In order to estimate molecular line emission, we have de- 
veloped a three dimensional non-local thermodynamic equi- 
librium (LTE) radiative transfer code based on an improved 
version of the Bernes (1979) algorithm. Our improvements 
(described more fully in Narayanan et al. 2006a,b) focus on 
including a subgrid model for giant molecular clouds (GMCs) 
in order to more accurately model the strongly density- 
dependent collisional excitation rates within our ~10 2 pc grid 
cells. We model GMCs as singular isothermal spheres (SISs) 
with power law index 2 (Walker, Adams & Lada, 1990), and 
assume half the cold gas mass in each cell is bound in GMCs. 
Our results for higher lying molecular levels and high dipole 
moment molecules are particularly improved using this sub- 
grid approach. 

The emergent spectrum is built by integrating the equation 
of radiative transfer through numerous lines of sight. The 
source functions in each grid cell are determined by the den- 
sities of molecules at levels u and / for a given transition 
u —*■ I. These level populations, «„,«/ are dependent on the 
incident radiation field from other clouds; we thus guess at a 
solution, calculate the mean intensity field in a Monte Carlo 
manner, determine the updated level populations through the 
rate equations, and repeat the process until the populations 
have converged. In between each iteration, the code goes into 



a sub-resolution process in which it decomposes each cloud 
into a SIS and determines the level populations in each sub- 
resolution element. The non-LTE aspect of this treatment is 
particularly important as the assumption of LTE breaks down 
when considering the propagation of radiation through me- 
dia with densities lower than the transition's critical density. 
For the radiative transfer calculations, we typically emitted 
~1 x 10 7 model photons per iteration, and consider the 2.73 K 
microwave background as the boundary condition. 

3. RESULTS 
3.1. Intensity Contour Maps 

In Figure [2 we show a series of snapshots in CO (J=l-0) 
intensity contour maps from the two merger simulations. The 
plot spans 45 h _1 Myr , and begins when the progenitors are 
approaching final coalescence. The black hole accretion rate 
nears its peak as the black holes merge at T^-4.15 h _I Gyr , 
and thus the feedback energy input from the AGN is near its 
maximum in the model with BHs. 

Beginning from Tr^I.13 h _I Gyr onward, the CO morphol- 
ogy of the galaxy in the BH model undergoes dramatic 
changes owing to feedback from the buried AGN. Massive 
blobs of cold molecular gas entrained in the wind are visi- 
ble through the CO (J=l-0) tracer. Indeed, while the AGN 
wind may in detail evaporate cold clouds via thermal conduc- 
tion, enough cold gas survives to be visible through molecular 
emission. Moreover, the clouds entrained in these outflows re- 
main cold and dense enough to continue forming stars. The 
dense cores in these clouds emit at CO transitions with rel- 
atively high critical densities; consequently, the outflows are 
visible against the background at transitions as high as CO 
(J=6-5). In order to produce the observed emission compara- 
ble to that of the nucleus, the outflows must have large col- 
umn densities. We find columns ranging from 5 x 10 22 cm~ 2 < 
7V(H 2 ) < lxl0 23 cm -2 through the outflow in Figure ^ de- 
pending on the viewing angle (although, we did not include a 
UV background in our models, and thus the true column may 
be less). As the outflowing gas leaves the nuclear region at 
velocities of ^200-300 km s" 1 , it becomes more diffuse, re- 
sulting in weaker CO emission. We find that the existence of 
molecular outflows is not unique to this particular model, and 
is seen in other merger simulations which include BHs. 

In the model without BHs (right column, Figure Q, the 
molecular outflows on the ^kpc scale observed in the BH 
model are not seen, highlighting the effects of AGN feedback 
in expelling loosely bound circumnuclear molecular gas. We 
note, however, that supernovae-driven winds are not incorpo- 
rated in these models. It is known, through absorption line 
spectroscopy, that such winds in starbursts can induce out- 
flows of comparable speeds (e.g. Heckman et al. 2000; Mar- 
tin, 2005; Rupke et al. 2005). While large columns of out- 
flowing molecular gas have not been imaged in many of these 
systems, there are notable exceptions such as the classic star- 
burst M82 (Walter, WeiB & Scoville, 2002). We will explore 
the effects of supernovae winds on the molecular gas in due 
course. 

3.2. Line Profiles 

The AGN induced outflows also leave their imprint on spec- 
tral line profiles. In Figure 13 we have calculated the spec- 
tral line emission from T= 1.155 h _1 Gyrin the model with 
black holes (lower left panel in Figure 0. The spectrum is 
generated along three orthogonal lines of sight. In order to 
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FIG. 1 . — CO (J=l-0) emission intensity contours for an equal mass galaxy 
merger. The left column is a time sequence for the model with BHs, and 
the right column without BHs. In the model with BHs, after the galaxies 
merge, massive amounts of gas are driven into the nuclear regions and accrete 
onto the supermassive black hole. The subsequent AGN feedback energy can 
blow large blobs of molecular gas out from the nuclear regions (see arrows). 
These features are not seen in the model without BHs. The time stamp of 
the image is in the top left of each panel, and is in units of fr'Gyr . The 
color contours are in units of K-km s (velocity-integrated Rayleigh- Jeans 
temperature), and the scale is on the right side of the plot. Each panel is 12 
hr'kpc on a side. The images are made at 1/4 h~'kpc spatial resolution. The 
coordinate system for use with Figure|2]is in the top right panel. 



simulate an unresolved observation, we set the merger at z=2 
(J7a=0.7, $7^=0.3, /z=0.75), and convolved the emission from 
the 12h _1 kpc image with a circular 30" (~235 kpc) Gaussian 
beam. 

Typically, in the BH model, once the galaxies have coa- 
lesced (when the BHs of the progenitors are indistinguishable 
in our simulations), the emission from the unresolved object 
is characteristic of a single Gaussian, centered at the systemic 
velocity of the galaxy. However, when viewing outflows with 
a strong line of sight (LOS) velocity component, a secondary 
peak appears superposed on the Gaussian emission line from 
the galaxy (e.g. left panel, Figure|2j; this peak is the emission 
from the outflow, and is redshifted or blueshifted from line 
center at the LOS velocity of the outflow. While we have pre- 
sented the CO (J=l-0) emission, these line profiles are similar 
through CO (J=6-5). 

The emission peak corresponding to the outflow appears at 
the greatest offset velocity with respect to the systemic veloc- 
ity of the galaxy when the outflow is moving mostly along the 
LOS (e.g. the -z observation, Figure 0. If the observation 
is tilted such that a smaller component of the outflow veloc- 
ity is along the LOS, the emission peak corresponding to the 
outflow will move closer to, and eventually will merge with 
the broad emission peak of the galaxy (y observation and -i 
observation, Figure|2] respectively). The peak temperature of 
the outflow emission may also decrease when the observation 
is tilted as the observer looks through less column. A given 
outflow along a particular line of sight is typically visible via 
its line profile for an average of ~10 h _1 Myr before the col- 
umn density through the outflow drops such that its emission 
is no longer detectable against the broader Gaussian emission 
from the central region. We estimate the "outflow" profiles 
are visible ~25% of the time in our simulations (which span 
200 h _I Myr ), averaged over many viewing angles. 

Double-peaked profiles have been observed in mergers 
which do not necessarily correspond to outflows: for example, 
observations of high density gas in the prototypical ULIRG, 
Arp 220, have evidenced a symmetric double-peak profile 
where each peak corresponds to the starburst regions of nu- 
clei of the progenitor galaxies (Taniguchi & Shioya, 1998, 
Sakamoto et al. 1999, Narayanan et al. 2005). Similar pro- 
files also occur in systems in which there is significant rota- 
tion, i.e. a disk galaxy or rotating nuclear ring. However, 
some differences exist between the double peaks originating 
in progenitor galaxies or rotating systems, and those caused 
by outflows. The double peaks characteristic of the former 
two cases are typically both broad, and symmetric about the 
systemic velocity of the galaxy. Conversely, the peak arising 
from the outflow is typically much narrower than the broad 
emission profile of the galaxy, due to its small velocity dis- 
persion along the LOS. That said, the component of the pro- 
file from the outflow can be quite bright, owing to large H2 
column densities through the outflowing material. 

The characteristic outflow profile consisting of a broad 
Gaussian with a narrow line superposed is also degenerate 
with that of high velocity gas falling in toward the nucleus. 
The infalling clouds tend to have a higher velocity dispersion 
than clouds entrained in the AGN wind by, on average, ~ 
30%, resulting in broader CO lines. It is unclear, though, 
from our simulations whether there is a significant enough 
difference between the velocity dispersion or column density 
in outflowing and infalling gas to determine the direction of 
flow from CO observations alone. There may be observational 
tests at other wavelengths to help break this degeneracy, how- 
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FIG. 2. — Synthetic CO (J=l-0) (rest-frame) emission line profile taken at T=1.155 fr'Gyr in the merger run with black holes over 3 orthogonal viewing angles. 
The direction of each line of sight is above the panels, and corresponds with the coordinate system in the top right panel of FigureQ The high velocity peak in 
the leftmost panel owes to the outflowing gas viewed with a large line of sight velocity component. In the y direction, the outflow has a smaller LOS velocity, 
and the emission owing to the outflow blends more into the main emission. In the -x viewing angle, outflows have a negligible effect on emission line. The best 
fit to the model spectra are overplotted with dashed lines. The spectra is taken such that the object is at z =2, and convolved with a 30" circular Gaussian beam. 



ever. In our simulations, infalling gas primarily leaves its im- 
print on the line profile just as the major merger has occurred, 
prior to the major AGN feedback phase. By serving as an 
indicator for AGN activity, X-ray observations may help dis- 
tinguish the origin of the "outflow" line profile. For example, 
using similar simulations to those explored in this study, Cox 
et al. (2006) have shown that the X-ray luminosity from dif- 
fuse gas in galaxy mergers peaks during the phase of heavy 
black hole accretion, and that thermal energy input from AGN 
feedback can produce X-rays consistent with observations of 
ULIRGs with known embedded AGN (e.g. UGC 5101, Iman- 
ishi et al. 2003). Elevated X-ray emission, both from diffuse 
gas, as well as hard X-rays from the central black holes may 
be indicative of a buried AGN. Correlations between hard X- 
ray flux and CO line profiles will be discussed in greater detail 
in Narayanan et al. (2006b). 

At z«2, where the infrared and submillimeter background 
are dominated by LIRGs and ULIRGs (Smail, Ivison & Blain, 
1997), and the quasar density is near its peak (e.g. Schneider 
et al. 2005), large molecular line surveys may prove fruit- 
ful in investigating the existence of a correspondence between 
double-peaked line profiles and quasar activity. Indeed, large 
fractions of sub-mm galaxies at high redshift show double- 
peaked line profiles (~ 50%, Greve et al. 2005, Tacconi et al. 
2006), some of which appear to have line profiles similar to 
those presented in the leftmost and middle panels of Figure|2j 
similarly, large fractions (~ 75%) of sub-mm galaxies show 
AGN activity as evidenced by IR and X-ray studies (Alexan- 



der et al. 2005, Polletta et al. 2006). 

4. SUMMARY AND CONCLUSIONS 

We have discussed two features of CO emission from cold 
molecular gas entrained in winds resulting from AGN feed- 
back in galaxy mergers: 1.) an extended CO morphology 
when the outflow is largely in the plane of the observation, 
and 2.) kinematic features in the emission line profile when 
the outflow has a significant line of sight velocity component. 
There may, of course, be hybrid cases as well, in which fea- 
tures of both signatures of AGN feedback are detectable. 

Observations of signatures such as these can help to inter- 
pret current observations of ULIRGs both at low and high-z. 
Emission maps and line profiles similar to those presented in 
this work may have already been observed in local mergers 
(e.g. NGC 985, Appleton et al. 2001), as well as in z « 2 sub- 
mm selected sources (Greve et al. 2005, Tacconi et al. 2006). 
Upcoming surveys with high resolution interferometers have 
the potential of constraining models of black hole growth and 
associated AGN feedback in galaxy mergers. 
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